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Abstract. We present a statistical analysis of the X-ray bursts observed from the 2002 June 1 8 outburst of the Anomalous 
X-ray Pulsar (AXP) IE 2259-1-586, observed with the Proportional Counter Array aboard the Rossi X-ray Timing Explorer. 
We show that the properties of these bursts are similar to those of Soft Gamma-Repeaters (SGRs). The similarities we find 
are: the burst durations follow a log-normal distribution which peaks at 99 ms, the differential burst fluence distribution is well 
described by a power law of index — 1 .7, the burst fluences are positively correlated with the burst durations, the distribution 
of waiting times is well described by a log-normal distribution of mean 47 s, and the bursts are generally asymmetric with 
faster rise than fall times. However, we find several quantitative differences between the AXP and SGR bursts. Specifically, 
there is a correlation of burst phase with pulsed intensity, the AXP bursts we observed exhibit a wider range of durations, the 
correlation between burst fluence and duration is flatter than for SGRs, the observed AXP bursts are on average less energetic 
than observed SGR bursts, and the more energetic AXP bursts have the hardest spectra - the opposite of what is seen for SGRs. 
We conclude that the bursts are sufficiently similar that AXPs and SGRs can be considered united as a source class yet there 
are some interesting differences that may help determine what physically differentiates the two closely related manifestations 
of neutron stars. 



INTRODUCTION 

Soft gamma repeaters (SGRs) are an exotic class of 
Galactic sources that are now commonly accepted as be- 
ing magnetars - isolated, young neutron stars that are 
powered by the decay of an ultra-high magnetic field. 
The evidence for high surface fields 10^'* — lO'^ G) 
comes from several independent lines of reasoning fHE 
|3l0]. These include: the high dipolar magnetic fields im- 
plied by the spin properties of SGRs seen in quiescence 
under the assumption of magnetic dipole braking 01^; 
the requirement of a magnetar-strength field to confine 
the energy released in the tails of hyper-Eddington out- 
bursts seen from two SGRs j|7| Isf]; the requirement of 
a high field to allow the decay rate necessary to power 
the burst and persistent emission 0,01; and the magnetic 
suppression of the Thomson cross-section, which allows 
hyper-Eddington bursts to be observed f?]. For a review 
of SGRs, see Kouveliotou et al. [these proceedings]. 

Anomalous X-ray pulsars (AXPs), another exotic 
class of Galactic neutron stars, have also been suggested 
to be magnetars |4|. This is because of their anoma- 
lously bright X-ray emission which can be explained 
neither by conventional binary accretion models nor ro- 



tation power fW\. Also, their spin parameters, as for 
SGRs, imply large magnetic fields under standard as- 
sumptions of magnetic braking. They also have similar, 
though on average softer. X-ray spectra compared with 
those of SGRs in quiescence. However, unlike SGRs, in 
the > 20 yr since the discovery of the first AXP |11|, 
none was seen to exhibit SGR-like bursts. For this rea- 
son, alternative models involving unconventional accre- 
tion scenarios have been proposed to explain AXP emis- 
sion See Kaspi et al. [these proceedings] for 
a review of AXPs. 

The magnetar model for AXPs was recently given 
a boost when SGR-like bursts were detected from two 
AXPs. Gavriil et al. lUsIl reported on the discovery of two 
X-ray bursts in observations obtained in the direction of 
AXP IE 1048.1-5937. The temporal and spectral prop- 
erties of those bursts were similar only to those seen only 
in SGRs. However, the AXP could not be definitely iden- 
tified as the burster. On 2002 June 18, a major outburst 
was detected unambiguously from AXP IE 2259H-586, 
involving over 80 bursts as well as significant spectral 
and timing changes in the persistent emission 1 16]. Those 
bursts demonstrated that AXPs are capable of exhibiting 
behavior observed, until now, uniquely in SGRs, there- 
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FIGURE 1. Distribution of the pulse phases of IE 2259+586 
which correspond to the times of the burst peaks (solid points). 
The solid curve is the folded 2-60 keV light curve of the 2002 
June 18 observation with the bursts omitted. 

fore implying a clear connection between the two source 
classes. Such a connection was predicted only by the 
magnetar model j3l • 

RESULTS 

The results presented here were obtained using the Pro- 
portional Counter Array [PCA; UtIi on board the Rossi 
X-ray Timing Explorer {RXTE). We use RXTE to moni- 
tor all five known AXPs on a regular basis as part of a 
long-term monitoring campaign [see 18, and references 
therein]. On 2002 June 18, during one of our regular 
monitoring observations, the AXP IE 2259H-586 exhib- 
ited an SGR-like outburst [iF]. The bursting behavior 
was detected by online RXTE monitors during the obser- 
vation. The observation spanned three orbits and had to- 
tal on-source integration time 10.7 ks. Following the out- 
burst. Target of Opportunity observations of the source 
were initiated the next day and continued at different 
intervals over the subsequent weeks, however no more 
bursts were seen. 

Our burst searching algorithm returned 80 significant 
bursts from the 2002 June 18 observation. Most bursts 
were single-peaked and had durations <1 s. A small 
handful (^12) were bright and had clear fast-rise, ex- 
ponential decay morphology. Four bursts were multi- 
peaked. Some bursts (^^5%) were approximately sym- 
metric, a few (^3%) fell faster than they rose while most 
fell slower than they rose. 

We calculated the occurance in pulse phase for each 
burst using the time of the burst peak and the rota- 
tional ephemeris given by Kaspi et al. |16]. Compar- 
ing the burst phase distribution to the pulse profile of 



FIGURE 2. Distribution of T90 durations for the bursts ob- 
served from IE 2259+586. The solid histogram line shows 
the observed binned distribution, while the dashed histogram 
line shows the corrected distribution (see text). The solid curve 
represents the best-fit log-normal model for the observed data, 
while the dashed curve is the best-fit log-normal model for the 
corrected data. 



IE 2259+586 at the time of the outburst, a correlation is 
seen (Fig.[0, where most of the bursts tend to occur when 
the pulsed intensity is high. We note that the two bursts 
seen from the AXP IE 1048.1-5937 |15| were also co- 
incident with the pulse peak, which strengthens the argu- 
ment that IE 1048.1-5937 was the source of those bursts. 

The T90 duration is the time between when 5% and 
95% of the total background-subtracted burst counts 
have been accumulated [e.g. il9.1 . SGR T90 distributions 
follow a log-normal distribution whose mean and stan- 
dard deviation vary with source [e.g."!?"]. At first we fit 
the measured values of T90 for the IE 2259+586 bursts 
with this model and found it to characterize the distribu- 
tion well (Fig. For low signal-to-noise bursts, T90 can 
be substantially underestimated. We corrected for this 
problem via Monte Carlo simulations [see 20). The cor- 
rected T90 distribution is shown in Figure|2] The best-fit 
mean is 99.31 ms with a range of 14.4-683.9 ms for one 
standard deviation. 

We measured the fluence of each burst and then 
grouped them in equispaced logarithmic bins. The dis- 
tribution of burst fluences is displayed in Figure |3] The 
low-end fluences are underrepresented because of sensi- 
tivity drop-off. Excluding the points having fluence <20 
PCA counts, the distribution is well modeled by a sim- 
ple power law. Using least-squares fitting we find a best- 
fit power-law index of —0.7 ±0.1, which corresponds 
to a differential spectrum dN/dF °^ /r-i-7±o.i 
ences in the 2-60 keV band range from ^5 x 10^^' to 
'--^ 7 X 10^'^ erg cm^^. These imply burst energies in the 
range 5 x 10^^ to ^ 7 x 10^'' erg, assuming isotropic 
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FIGURE 3. Distribution of the 2-60 keV fluence for each 
burst. Solid points represent average values of fluence in eq- 
uispaced logarithmic bins for which our observations had full 
sensitivity. The open points suffered from reduced sensitivity. 
The best-fit line was determined using the solid points only and 
is shown as a solid line; the dashed lines are its extrapolation. 

emission and a distance of 3 kpc to the source f2l"].The 
sum total of all burst fluences is 5.6 x 10^^ erg cm^^, 
corresponding to energy 6.0 x 10^' erg (2-60 keV). 

Gogu§ et al. 1 19] also find a clear correlation between 
burst durations and total burst fluence. A similar correla- 
tion is seen in our IE 2259H-586 data. To quantify it, we 
grouped the T90 values in equispaced logarithmic bins 
and determined group-averaged fluences for each bin. 
Least-squares fitting to a simple power-law model yields 
F oc T+5'■^'^±''■''^ with reduced = 1-0. 

We measured the peak flux in a 61.25-ms time bin 
for each burst. Our burst-identifying algorithm is less 
sensitive to bursts of smaller peak flux, we compen- 
sated for this effect via simulations [see 20]. Using 
least-squares fitting we found that the corrected distri- 
bution is well modelled by a simple power law with in- 
dex — 1.42 ±0.13. Peak fluxes in a 61.25-ms time bin 
range from ^ 1 x 10^^ to ^ 1 x 10^^ erg cm^^ s^', 
which imply peak luminosities in the range ^ 1 x 10^^ 
to ^ 1 x 10^^ erg s^'. On shorter time scales we find 5 
bursts with peak fluxes which are super-Eddington. The 
peak fluxes in a 1/2048 s time bin for these bursts range 
from ~ 2 X 10^^** to - 8 x 10^^ erg s^^ 

Burst rise and fall time distributions were well mod- 
elled by log-normal distributions. For the rise time distri- 
bution, we find a mean of 2.43 ms and a range of 0.51- 
11.51 ms for one standard deviation. For the fall time 
distribution, we find mean 13.21 ms and a range of 3.52- 
49.55 ms for one standard deviation. 

SGR waiting times, defined as the temporal separa- 
tions of adjacent bursts, are found to follow log-normal 
distributions i2^E3ll . We measured the waiting time for 



FIGURE 4. Hardness ratio (H) versus fluence (F). The open 
points are hardness ratio measurements for individual bursts. 
The solid points are weighted averages of hardness ratios for 
bursts in equispaced logarithmic fluence bins. The line repre- 
sents the best-fit logarithmic function for the weighted aver- 
ages, /f = 0.31 x log F - 0.09. 



the IE 2259H-586 events, excluding those interrupted by 
Earth occultations. Our waiting time distribution is also 
well modelled by a log-normal distribution with mean 
46.7 s and a range of 10.5-208.4 s for one standard de- 
viation. We find no correlation between the burst energy, 
duration and the waiting time until the next burst, nor 
with the elapsed time since the previous burst. There is a 
correlation between the waiting time and the event time, 
which implies that the mean of our waiting time distribu- 
tion depends on the time at which we started observing 
the outburst. 

Gogu§ et al. fl^ noted that SGR bursts tend to soften 
with increasing burst energy. We studied the hardness ra- 
tio/fluence relationship by extracting spectra and creating 
response matrices separately for each burst. Hardness ra- 
tios were defined as the ratio of the counts in the 10- 
60 keV band to those in the 2-10 keV band as in Gogu§ 
et al. 1 19]. Also following Gogiig et al. ] 19], we divided 
the bursts into equispaced logarithmic fluence bins and 
calculated a weighted average hardness ratio for each 
bin. Figure |3 shows the weighted mean hardness ratios 
as a function of fluence. A clear positive correlation is 
seen. 



DISCUSSION 

The bursts we have observed for IE 2259H-586 are clearly 
similar to those seen uniquely in SGRs. As concluded by 
Gavriil et al. J 15] and Kaspi et al. fld\, AXPs and SGRs 
clearly share a common nature, as has been predicted by 
the magnetar model. The bursts seen in the 2002 June 



18 outburst of IE 2259+586 are qualitatively similar to 
those seen in SGRs, and in many ways quantitatively 
similar. Specifically: 

• the mean burst durations are similar 

• the differential burst fluence spectrum is well de- 
scribed by a power law of index —1.7, similar to 
those seen in SGRs (and earthquakes and solar 
flares) 

• burst fluences are positively correlated with burst 
durations 

• the distribution of and mean waiting times are sim- 
ilar 

• the burst morphologies are generally asymmetric, 
with rise times usually shorter than burst durations 

However, there are some interesting quantitative differ- 
ences between the properties of the AXP and SGR bursts. 
These may help shed light on the physical difference(s) 
between these classes. The differences can be summa- 
rized as: 

• there is a significant correlation of burst phase with 
pulsed intensity, unlike in SGRs 

• the AXP bursts have a wider range of burst duration 
(though this may be partly due to different analyses 
procedures) 

• the correlation of burst fluence with duration is flat- 
ter for AXPs than it is for SGRs (although when 
selection effects are considered, this correlation 
should really be seen as an upper envelope for AXPs 
and SGRs) 

• the fluences for the AXP bursts are generally 
smaller than are in observed SGR bursts 

• the more energetic AXP bursts have the hardest 
spectra, whereas for SGR bursts, they have the soft- 
est spectra 

• under reasonable assumptions, SGRs undergo out- 
bursts much more frequently than do AXPs 

Given the rarity of AXP bursts coupled with the 
unique information that detection of such bursts pro- 
vides, observing more outbursts is obviously desirable. 
Continued monitoring is thus clearly warranted, and 
RXTE with its large area and flexible scheduling is the 
obvious instrument of choice. 
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